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ABSTRACT 

We present a model that describes stellar infrared excesses due to heating of the interstellar (IS) 
dust by a hot star passing through a diffuse IS cloud. This model is applied to six A Bootis stars with 
infrared excesses. Plausible values for the IS medium (ISM) density and relative velocity between the 
cloud and the star yield fits to the excess emission. This result is consistent with the diffusion/accretion 
hypothesis that A Bootis stars (A- to F-type stars with large underabundances of Fe-peak elements) 
owe their characteristics to interactions with the ISM. This proposal invokes radiation pressure from 
the star to repel the IS dust and excavate a paraboloidal dust cavity in the IS cloud, while the metal- 
poor gas is accreted onto the stellar photosphere. However, the measurements of the infrared excesses 
can also be fit by planetary debris disk models. A more detailed consideration of the conditions to 
produce A Bootis characteristics indicates that the majority of infrared-excess stars within the Local 
Bubble probably have debris disks. Nevertheless, more distant stars may often have excesses due to 
heating of interstellar material such as in our model. 

Subject headings: stars: A Bootis stars - circumstellar matter - stars: infrared excesses 



1. INTRODUCTION 

The vast majority of stars follow simple trends in mass, 
luminosity, age, metallicity, and other parameters as cap- 
tured in the simplest form through the Russel-Vogt The- 
orem. The exceptions to these trends therefore must 
represent interesting events in stellar evolution. A Boo- 
tis stars are an outstanding example. They are late B- 
to early F-type Population I stars with peculiar surface 
abundances: while light elements like C, N, O and S are 
roughly solar, heav y elements (Fe-peak) are depleted up 
to a factor of 100 (jPaunzen et al.ll2002f ). The ages of 
A Bootis stars span the entire main sequence lifetimes 
for their types, with a distribution that has a peak at 
a rather evolved stage (~ 1 Gyr). Less than 2% of all 
objects in the relevant spectral domain are believed to 
belong to the A Bootis group, suggesting either a very 
short time scale (10^ years) for the phenomenon, or un- 
common conditions for its occurrence. 

Understanding A Bootis stars will provide important 
insights to more general problems. Due to their ex- 
tremely shallow convection zones, A-type stars in gen- 
eral and chemically peculiar ones in particular provide 
excellent natural laboratories for the study of fundamen- 
tal astrophysical processes such as diffusion, meridional 
circulation, stellar winds, and accretion of circumstellar 
a nd interstellar materia l. 

iMichaud k, Charlandl ()1986l ) suggested that the pecu- 
liar chemical abundances on the surfaces of A Bootis 
stars are due to accretion of circumstellar material that 
is mixed in the shallow convection zone of the star by the 
joint action of gravitational settling and radiative accel- 
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eration. ICharbonneaul (1X9931) [Venn fc LambertI (|1990( l. 
and iTurcotte &: Charbonneaul (|1993f ) developed this hy- 
pothesis further. It naturally explains why the anoma- 
lous abundance pattern is similar to that found in the 
gas-phase of the interstellar medium (ISM), where metal- 
lic elem ents like iron and silicon h ave co ndensed into dust 
grains. lAndrievskv fc PaunzenI ()2000f ) suggested a spe- 
cific mechanism: that A Bootis behavior results from 
condensation of dust grains in the circumstellar envi- 
ronment, their removal through radiation pressure and 
the accretion of the remaining metal-poor gas. However, 
their diffusion/ accretion model has many free parame- 
ters and fails to explain some observations such as the 
small proportion of stars with A Bootis characteristics. 
Alternative theories in voke contact binar y systems merg- 
ing into a single star (lAndrievskvlll997|) and unres olved 
spectroscopic binarie s IParaggiana fc Bonifaciolll999 ) . 

More recently, iKamp fc PaunzenI ( 20021 ) have 
suggested that the A Bootis phenomenon re- 
sults from interactioii along the lines discussed by 
lAndrievskv fc PaunzenI (2000) but between the star and 
the local ISM (Figure [1]). This new hypothesis removes 
the ad-hoc assumptions and some of the degeneracy of 
previous models. Different levels of under abundance 
are produced by different amounts of accreted material 
relative to the photospheric mass. The small fraction of 
A Bootis stars is explained by the low probability of a 
star-cloud interaction and by the effects of meridional 
circulation, which washes out any accretion pattern 
a few million years aft er the accretion has stopped 
(|Kamp fc PaunzenI I2OO2D . The behavior is suppressed 
in low temperature stars because of their more massive, 
difficult-to-contaminate convection zones. Strong stellar 
winds prevent the accretion of material for very hot 
stars. 

In addition to their peculiar abundance pattern, some 
A Bootis stars show infrared excesses, suggestive of the 
presence of dust. Since any ongoing interaction with the 
ISM will heat interstellar dust, these cases can probe 
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Fig. 1. — Geometry of the star-cloud interaction. Two differ- 
ent avoidance radii, corresponding to two different grain sizes are 
shown. 



the model propos ed by Kamp fc Paunzenl ()2002D . Using 
IRAS photometry. iKin^ ( 19941 ) found that 5 out of a sam- 
ple of 56 A Bootis stars have infrared ex cess in at least 
one of the satellite's bands. More recently. iPaunzen et al.l 
(|2003f ) have compiled data from the ISO and IRAS satel- 
lites to identify A Bootis stars with excess infrared emis- 
sion. According to their work, 23% of 26 well established 
members of this group show infrared excesses with de- 
rived fractional dust luminosities (Lir/L*) ranging from 
2.2 X 10~^ to 4.3 X 10~*. Both the incidence of excesses 
and the fractional luminosities are similar to the char- 
acteristics of planetary debris disks for stars of similar 
s pectral type (iRieke et al.l l2005l: ISu et al.ll2006l ) . 

iKamp et al.l ( 2002D us"ed a sample of 21 A- type stars 
including debris-disks, dust-free stars and six A Bootis 
stars, and found no fundamental difference between the 
atmospheric compositions of the stars with and without 
debris disks. This result suggests that circumstellar ma- 
terial is unlikely to be linked to the anomalous A Bootis 
abundances. Thus, the IR excesses of the A Bootis-type 
stars might be due to a different mechanism, e.g., ac- 
cretion of interstellar material distributed in a diffuse 
cloud. This i dea is also support ed by coronagraphic ob- 
servations by Ik alas et al.l (|2002f ) that showed some stars 
previously catalogued as debris disk stars are surrounded 
by an extended reflection nebulosity. 

In this paper we develop a the rmal emission model 
to explore the iKamp fc Paunzenl ([2002) ISM interac- 
tion hypothesis. We emphasize the basic physics behind 
the model rather than detailed simulations of individual 
stars. Our most important conclusion is that simple in- 
frared excess SEDs can often be fitted equally well by 
either a debris disk or interstellar interaction model. As 
a result, further probing of ISM-related excesses requires 
detailed modelling of stars where other lines of evidence 
positively identify such interactions (e.g., for S Velorum, 
iGaspar et a"ni2008f ). 

In §2 we provide a description of the observations used 
in this paper and discuss the physical parameters of our 
sample of stars. We introduce our thermal emission 
model in §3, where we describe the basic physical param- 
eters involved in the calculation of the spectral energy 
distributions. §4 explores the effect of varying different 
model parameters on the output SEDs. We discuss the 
implications of our results for the identification of debris 
disks and the nature of A Bootis stars in §5 and conclude 
the paper in §6. 



2. SAMPLE SELECTION AND OBSERVATIONS 

Our sample consists of six A-ty pe stars with well es- 
tablished A-Bootis characteristics (jPaunzen "erall[200l) 
and at distances less than 100 pc. We utilize pho- 
tometry at 24 and 70 of these stars obtained with 
the Multiband Im aging Photometer for Spitzer (MIPS; 
iRieke et al.l2004l) and reduce d with the Instrumen t Team 
Data Analysis Tool (DAT; iGordon et 311120051 ). The 
stars were required to have 24 /xm excesses at least 10% 
above the photospheres and detected excesses also at 70 
/im. We have excluded binary stars (e.g., HD 198160), 
where the companion may contaminate the photometry. 

Details of the reduction proced ures and calibration are 
provid ed in iGordon et al.l ()2007f ) and lEngelbracht et all 
()2007l ). Flux densities were extracted using aperture 
photometry with radii of 14.94" at 24 /im and 35" at 
70 nm, relative to sky annuh at 29.88" to 42.33" and 39" 
to 65" , respectively. Aperture corrections appropriate to 
a point source were applied to the photometry. The pho- 
tometric aperture was selected by consider ation of the re- 
solvcd structure at 24 and 70 fim for 6 Vol (jGaspar et al.l 
[2008), where the bulk of the detected flux would fall 
within our selected apertures. Since our six stars are all 
at greater distances than 6 Vel, we should have captured 
virtually all the flux detectable with MIPS. 

To determine the excess properties, we subtracted the 
photospheric fluxes from the measurements. We de- 
termined these fluxes using Kurucz atmospheric models 
matched to the effective temperature and surface grav- 
ity for each star (assumed to have solar abundance). 
Stellar masses were obtained f rom the literature (e.g. 
lAllende Prieto fc Lambertlll999[ ) and they are represen- 
tative of the spectral type of each star. Stellar radii 
are fine-tuned to produce a photospheric SED that is in 
agreement with optical and near-infrared (2MASS) pho- 
tometry, while total luminosities were calculated by in- 
tegrating the resulting energy distribution of the scaled 
Kurucz models over the wavelength range where they are 
defined. Table 1 lists the fundamental stellar parameters 
as well as the predicted and measured fluxes at 24 and 
70 /im. 

3. MODEL DESCRIPTION 



Following the iKamp fc Paunzenl (2002^ interaction 
model, here we consider an A-type star moving through 
a uniform density, optically thin ISM cloud. The relative 
motion between the star and the cloud is described by a 
constant relative speed v^ci- The emission of the central 
star irradiates the sub-micron dust grains present in the 
cloud, warming them up to equilibrium temperatures, at 
which they re-emit the energy in the IR. Since the ISM 
cloud is optically thin, thermal radiation emitted by the 
warm grains will not be reprocessed by the cloud itself. 
Our model calculates the resulting SED at mid-infrared 
wavelengths. 

3.1. Geometry of the cloud 

Our model assumes a spherical cloud, penetrated by an 
axially symmetric paraboloidal cavity excavated by the 
radiation pressure of t he star as it moves thro u gh th e 
cloud, as described by lArtvmowicz fc ClampinI (|1997f) . 
We can describe this geometry mathematically using 
spherical coordinates (r, 6, 4>). Because of the axisym- 
metrical character of the cavity, only a radial coordinate 
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TABLE 1 

Physical parameters of the sample stars. 



HD No. 


d 


Toff 


logio L* 


R* 


M* 


F24 


F70 


Fpr24 


FprTO 


T24-70 




pc 


K 


Lq 


Rq 


M© 


mJy 


mJy 


mJy 


mJy 


K 


11413 


74.8 


7925(124) 


1.32 


2.4 


2.6 


53(1) 


56(2) 


47.71 


5.42 


80 


30422 


57.5 


7865(108) 


0.97 


1.6 


1.8 


47(2) 


65(2) 


36.21 


4.10 


76 


31295 


37.0 


8920(177) 


1.23 


1.8 


2.1 


177(3) 


438(28) 


101.66 


11.55 


71 


110411 


36.9 


8930(206) 


1.15 


1.6 


2.0 


147(4) 


256(10) 


91.66 


10.45 


84 


125162 


29.8 


8720(156) 


1.20 


1.9 


2.0 


286(3) 


392(6) 


190.74 


21.59 


86 


183324 


59.0 


8950(204) 


1.16 


1.8 


2.1 


53(2) 


30(5) 


43.70 


5.04 


84 



^ MIPS measurements and predicted photospheric fluxes are shown. Errors shown corr espond to 
measu rement uncertainties. Efl'cctive temperatures (including err o rs) we re taken from Pau nzen et al.l 
(HQOI), except HD 97633, taken from lAUende Prieto fc Lamberti | |T999|) . T24_7o refers to the 24-70 



jjLm color temperature. 

(r) and an azimuthal coordinate {0) are needed to de- 
scribe its surface, which is symmetrical with respect to 
4). We choose the coordinate system so that the star is at 
the origin (r = 0) and marks the center of the spherical 
structure. The ray with zero azimuthal angle coincides 
with the direction of relative motion (Figure [T]). The 
equation describing the surface of the cavity is then 



(1) 



1 + cos 6* 



The geometry is parametrised by the avoidance radius, 
Tin, which determines the shape and size of the cavity^. 
The grains never approach the star closer than rin. The 
avoidance radius can be estimated by considering energy 
conservation. For a relative velocity between the star and 
the cloud (t'rei)j a ratio of radiation pressure to gravity 
, and a stellar mass of M, , it can be show that 



2(/3- 1)GM* 



(2) 



where G is the Newtonian gravitational constant 
(lArtvmowicz fc ClamDinlfl997l) . 

iBurns et al.l (|1979l ) show how to calculate (3, which de- 
termines Tin as a function of the grain size a. We briefly 
reproduce their approach. The gravitational attraction 
of the star upon a spherical particle of radius a and den- 
sity pdust (g cm~^) that is located at a distance r from 
the star is 



= -7raVdustGM,/r2 



(3) 



The pressure force due to the radiation field of the star 
is given in terms of the stellar luminosity and the 
geometric cross section A of a spherical grain of radius 
a: 



(4) 



where is the average radiation pressure coefficient 
for grains of size a. To calculate this coefficient, we first 
compute the grain-size-dependant radiation pressure ef- 
ficiency as a function of wavelength: 

Qpr(A)=Q:bs + Qsca(l-(cOSa'^)). (5) 

^ The behavior is analogous to Rutherford scattering of posi- 
tively charged particles by a positively charged atomic nucleus. 



Both the scattering coefficient Qg^a ^^'^ ^^'^ absorption 
coefficient (5°^,s, as well as the scattering angle, a°, are 
functions of the wavelength. They are available for ap- 
propriate grain parameters in tabulated form. Once we 
have determined the shape of Qp^{X), we compute Qp^ as 
a weighted average of Q^^. over the wavelength spectrum: 



/FK(r.)Qg,.(A)dA 

/FK(r*)dA 



(6) 



where Fk(T,) is the energy density of the respective Ku- 
rucz model. We can then compare the radiation force to 
the gravitational force by means of their ratio: 



Z?-^ = F,/Fg = (3L,/16^GM,c)(Q^,/apdust) 

-1 / \ -1 



Pdust 

g cm~' 



(7) 



where and M, are in solar units, pdust should be dis- 
tinguished from the total dust mass density of the cloud, 
Pcioud, which is one of the free parameters in our model. 
Eq. (7) shows the dependence of /3° on the grain size. 
As a consequence, under the same physical conditions, 
grains of the same composition but different sizes will 
have different avoidance radii, which in turns leads to 
different equilibrium temperatures. Mathematically, this 
dependence arises not only from the factor in Eq. (7), 
but also because the optical properties of a grain, (i.e., 
Qahs ^'^^ Qsca) ^^'^ fuuctious of the grain size. In §4 
we will discuss further variations in rin due to different 
compositions of the grain material. Because of different 
grain properties, the model cavity has a large range of 
avoidance radii, resulting in a corresponding range in the 
cavity size and in the grain equilibrium temperatures. 

3.2. Relative velocities 

As shown in Eq. (2), the shape of the cavity depends 
on the velocity of the star relative to the cloud. De- 
termination of the heliocentric velocity for the interstel- 
lar material close to our stars, wism, and hence, of the 
relative velocity of the system {vre\), is not straightfor- 
ward. In some instances, there are candidate clouds with 
known velocity vectors in the local ISM which could be 
responsible for interacting with the A Bootis stars. A 
catalogue of local nearby clouds of warm ISM and their 
velocity vectors an d approximate outlines is presented in 
iRedfield fc Linskvl (f2008iV A program called "windsock" 
was written by Eric Mamajek which estimates the rel- 
ative 3D velocities of ISM clouds and the target stars. 
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and calculates the projected relative tangential and ra- 
dial motions appropriate for comparing to the detailed 
dust interaction models. The clo ud vel ocity vectors were 
adopted from Redfield fc Linskvl (|2008f ) , and the velocity 
vectors of the stars were estimate d from revised Hippar- 
cos astrometry from Ivan Leeuwe n (2007) and the c om- 
piled radial velocity catalogue of Gontcharov (2006') us- 
ing the formulae of Johnson & Soderblom (1987). Table 
2 shows the velocities for each of our stars relative to 
several candidate clouds, as calculated using windsock. 
Since we cannot determine in each case which of these 
clouds actually interacts with the star, we have decided 
to study the behavior of the model at two limiting veloc- 
ities, which cover approximately the range of velocities 
in Table 2. These velocities are 15 km s~'^ and 40 km 
s^^. It is important to mention that the relative veloc- 
ity for two of our stars (HD 125162 and HD 30422) is 
closer to our lower velocity limit (15 km s~^), while for 
the remaining five stars this velocity is closer to the up- 
per limit (40 km s""'^). In §4 we study the effect of the 
relative velocity on the infrared excess properties and we 
estimate the corresponding uncertainties. 



3.3. Thermal 



emission 



With the geometric distribution of the dust defined, we 
can predict the infrared fluxes in the MIPS bands. The 
fixed parameters in our model are the relative velocity 
Urei, the mass to luminosity ratio (M/L) of the star, its 
radius and heliocentric distance d, as well as the com- 
position, size and optical properties of the dust grains. 
All except the last set of parameters have already been 
discussed. 

For the grains, we assume an interstellar composi- 
tion, i.e. the sirn plistic homogeneous sphere models 
(jMathis et al.lll9770 . hereafter MRN, with two (graphite, 
silicate) components and a MRN size distribution of the 
form dn = Ca~^'^da with lower and upper limits at 
amin = 0.005 /im and amax — 0.25 /im respectively 
([Mathisl |199tf) . Graphite and silicate populations are 
treated independently. The free parameters in our model 
are then just the outer radius of the cloud rout, i-e., its 
size, and the total mass density of its dust component, 

Pcloud ■ 

To calculate IR fluxes, we assume that the grains are 
in radiative equilibrium. Grains absorb stellar radiation 
and re-emit it at temperature Tg. We use the Kurucz 
model described in §2 for the stellar radiation field. The 
radiative equilibrium at a distance r from the star can 
be expressed as 

J QlsWFK{X)dX = J Q^,,(A)B(A,rg)dA. (8) 

QahsW the absorption efficiency for grains of size a. 
Since we assume thermal equilibrium, Kirchoff's law ap- 
plies and the absorptivity of the dust grains equals their 
emissivity. Qabs thus appears on both sides of Eq. (9), 
but evaluated over different wavelength regimes. By solv- 
ing this equation numerically, we find the equilibrium 
temperature Tg as a function of distance, grain size and 
grain composition. 

The functional forms of (5abs('^) interstellar dust 
Kr a,ins of the des i red co mpo sition have been publ ished 
bv iDraine fc Led (|1984D and iLaor fc Draind (|1993f ). for 



grain sizes covering the range 0.001 /im < a < 10 /im 
and for wavelengths between 0.001 /im and 1000 /im. 
We use these tables in our integration of Eq. (8). For 
each grain size, a different avoidance radius is calculated 
according to Eq. (2). The grain number density of the 
cloud, n, is related to the grain size by 



n{a) = Ca ''•^da. 



(9) 



The integration is performed over the range of grain 
sizes. The astronomical silicates (pdust = 3.3 g cm~^) 
and graphite (pdust = 2.25 g cm~'^) both have a size dis- 
tribution described by Eq. (9), and each of them carries 
half the total mass (Mdust) of the dust component of the 
ISM cloud. Mathematically, this leads to the normalisa- 
tion equation: 



n{a)m{a)V {a)da 



M. 



dust 



dust 



(10) 



where n{a) and m{a) = An Pdust/3 are respectively 
the number density and individual grain mass as func- 
tions of the grain size, V{a) is the volume occupied by 
grains of size o (each grain size has a different avoidance 
radius) and Vdust is the total volume occupied by the 
dust cloud. Combining Eq. (9) and Eq. (10) and solving 
for the normalisation constant C, we get: 



C 



3/^ Vdust 

Pdust 



y(a)a-"-Ma 



(11) 



which is valid for any pdust- Taking into account Eq.(l) 
for each grain size, we have calculated the normalisation 
constants for each population of grains as well as the 
total dust mass of the cloud, Mdust- 

Now that we have characterised the grain size distri- 
bution, we can calculate the infrared flux. For each grain 
size, we integrate the contributions of the modified black- 
bodies (dust grains), as determined by the distribution 
of equilibrium temperatures (Eq.(8)) over the cloud vol- 



dFa(A) - n(a)^Q:b.(A) B(A, Tg(r))dV (12) 

Finally, we integrate over the grain size distribution to 
obtain the spectral energy distribution (SED): 



F(A) = 



i?(A,Tg(r))dV (13) 



We fit the MIPS photometry with this function. 

If the IR excesses of the A Bootis stars are indeed due to 
a thermally emitting nebulosity, the MIPS observations 
discussed in §2 may partially resolve the sources. In that 
case, the fitting to the photometry needs to allow for the 
source extent. We discuss this issue in §4.1. We then 
continue to describe the underlying model physics and 
estimates of cloud dimensions and density to assess the 
plausibility of the star-ISM interaction model. 

4. EXPLORING THE PARAMETER SPACE 

To assess the uncertainties involved in our model and 
to test its validity, we perform a series of trial calculations 
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TABLE 2 

WiNDSOCK MODEL VELOCITIES. 



HD No. -ULIC '"Vol fBIuc I'Aur fGom fNGP "Leo I'Mic ^^Aql ^Eri 

11413 40.0 43.4 58.0 - - - _ _ _ _ _ 

30422 12.2 - - 7.3 9.6 - _ _ _ _ _ 

31295 27.5 - - - 18.6 29.0 _ _ _ _ _ 

110411 _ _ _ - _ _ 56.8 43.3 

125162 _ _ _ _ _ - 2.3 10.4 

183324 - 42.0 _ _ _ _ _ _ 42.1 64.7 36.8 



For each star, velocities in km s ^ are shown for c louds that are more likely to be 
interacting with it according to lRedfield fc Linskvl 1)20081 ) . 



aimed at establishing the impact of individual parame- 
ters on the SED and derived quantities. We also examine 
the influence of grain composition and total dust density. 
For these tests, we use the stellar and cloud parameters 
derived for HD 125162, the prototype of the A Bootis 
stars. 

4.1. Outer radius of the cloud 

Because of its relatively small heliocentric distance, HD 
125162 is the most demanding example to assess the un- 
certainties due to our selected aperture for photometry. 
At the distance to this star (29.8 pc), the aperture radius 
of 14.94" used to perform the point source photometry 
at 24 /im corresponds to a physical radius of 445 AU, 
whereas an aperture radius of 35" , used to perform the 
photometry at 70 /Ltm, corresponds to a physical radius 
of 1043 AU. As discussed in the next section, to produce 
the observed fluxes, our model suggests an outer radius 
of 2200 AU for the diffuse interstellar cloud, when a ve- 
locity of 40 km s^^ is used. To estimate the uncertainties 
in our best fit parameters due to the difference between 
the aperture radius and the actual size of the cloud, we 
compare the fluxes at 24 and 70 /xm obtained for the 
entire cloud (rout = 2200 AU) with the corresponding 
fluxes that now only include material within the corre- 
sponding aperture radii (445 and 1043 AU). At 24 /im 
this flux is calculated in the same way as we obtained 
the photometric measurements, i.e., we subtract the flux 
from the sky annulus and apply an aperture correction 
factor: 

F(24 ^xm) = 1.143(i^2T"'*'"'° " i^aT"'"') (14) 

We obtain a flux that is only 20% smaller than the flux 
at 24 /im calculated using the entire cloud. At 70 /im, 
the emission from the sky annulus is comparable to that 
from the aperture and hence we did not subtract it. We 
compare the flux obtained for the entire cloud with that 
calculated with an outer radius of 1043 AU, correspond- 
ing to the 35" aperture. In this case, the flux turns out 
to be a factor of 2.5 smaller than the one obtained for 
the entire cloud. This result suggests that a larger pho- 
tometric aperture should have been employed. However, 
there are significant difficulties in identifying and mea- 
suring low surface brightness extensions at 70 /<m and it 
is not clear that we would improve the overall accuracy 
of the measurement. Since HD 125162 is the closest star 
in our target list, we consider these uncertainties to be 
an upper limit for our sample. 

Figure [2] illustrates a more general case. It shows dust 
SEDs for five diffuse clouds with different outer radii. To 



exclude other possible effects, we have used here only one 
composition (silicates) and a single grain size (0.008 /im, 
corresponding to the weighted mean of the distribution). 
The stellar parameters correspond to HD 125162, with a 
velocity corresponding to our higher limit of 40 km s~^. 
Even though, as we stated before, for HD 125162 the rel- 
ative velocity is closer to our lower limit of 15 km s~^, 
here we are only interested in the impact that varying 
the outer radius has on the model and hence the velocity 
used is not relevant. Additionally, as showed in §4.3, the 
SED does not change significantly in the mid-infrared re- 
gion within a range of velocities between 15 and 50 km 
s~^. The total dust mass density is kept constant, and 
the outer radius ranges from an upper limit of the typ- 
ical radii of circumstellar disks (~ 500 AU) up to the 
size scale of diffuse clouds in the interstellar medium 
10^ AU). As the cloud becomes larger, the total emis- 
sion increases and the maximum of the SED is shifted to 
longer wavelengths. As we increase the outer radius of 
the cloud we increase the contribution to the SED of ex- 
ternal, colder material, which dominates the total mass. 
Since the density is kept fixed, increasing values of rout 
imply an increase in the total cloud mass and therefore 
an increase of the total emission. A tend ency for ISM- 
interaction infrared excesses to be " warm" (jGaspar et al.l 
[200l might be exaggerated observationally by the in- 
herent measurement bias against this dilute, extended 
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Fig. 2. — Dependence of the dust IR emission on the outer radius 
of the cloud. The size range goes from typical disks sizes (500 AU, 
solid line) up to typical sizes of diffuse interstellar clouds (100000 
AU, triple dot-dashed line), with intermediate values at 3000 AU 
(dotted line), 10000 AU (dashed line) and 50000 AU (dot-dashed 
line). 
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Fig. 3. — Dependence of the IR emission on the relative star- 
cloud velocity. Velocities are: 0.9 km s~^ (solid line), 1.5 km s~^ 
(dotted line), 4 km s~^ (dashed hne), 10 km s~^ (dot-dashed line), 
35 km s~^ (triple dot-dashed line) and 100 km s~^ (long-dashed 
line) 

emission component. 

4.2. Relative velocity 

Given the stellar space velocities and the available mea- 
surements on the local interstellar cloud kinematics, the 
parameter v^^i can lie approximately between 15 km s~^ 
and 40 km s~^, as stated. Figure [3] shows the calculated 
SEDs for a cloud of silicate grains with single grain size 
a = 0.008 fj.m around HD 125162. The relative veloc- 
ity spans a range between 0.9 km s~^ and 100 km s~^, 
while all the other parameters remain fixed. Changing 
the relative velocity changes the avoidance radius accord- 
ing to Eq. (2). However, for the range of velocities of 
interest (15 km s~^ to 40 km s~^), the SEDs are very 
similar at wavelengths larger than 20 ^m, and approach 
an asymptotic SED. This asymptotic case corresponds 
to high velocities (wrci > 10 km s~^ for silicate grains), 
in which case rin = and the cavity vanishes, since there 
is no time to excavate it. It is good to remind ourselves, 
however, that physically this will not be the case, since 
there will always be an inner radius set by the subli- 
mation temperature o f the grains (~ 1500 K for silicate 
grains, 'Salpctcr 1977^. We see in Figure [3] again a color 
change of the SED similar to the one described in §4.1: 
as the relative velocity becomes smaller, the avoidance 
radius becomes larger and colder material dominates the 
emission. Carbonaceous grains behave similarly, but due 
to their larger radiation pressure efficiency the velocities 
needed to reach the limiting case are of the order of f rci > 
100 km s^^. As shown in Tables 3 and 4, for a particular 
star-cloud system, derived cloud densities at 15 km 
are approximately 3.5 larger than in the 40 km s~^ case, 
while derived outer radii are about 1.6 times smaller. 
The total mass needed to produce the measured excess, 
however, remains almost the same in both cases. 

4.3. Dust composition 

For the purposes of this paper, we use a mixed popu- 
lation of astronomical silicate and graphite dust grains, 
each of them providing 50% of the total dust mass. We 
assume bare grains, without ice mantles. The com- 



FlG. 4. — Dependence of the avoidance radius on the composition 
of dust grains. The two different compositions are: astronomical 
silicate grains (dotted line), and graphite grains (solid line). The 
encounter velocity is 15 km s~^, and the stellar parameters are 
those for HD 125162. 

plex refractive indices and the abs orption and scatter- 
ing properties of du st grains (e.g. iDraine &: Leel[l98^ 
iLaor fc Draind '1993) imply that at the wavelengths at 
which A-type stars emit most of their flux, carbona- 
ceous compounds are more efficient at absorbing UV 
photons compared to astronomical silicates, according 
to Eq. (6). Hence, the radiation pressure on graphite 
is larger than for the silicates and, as a consequence, car- 
bonaceous grains will be located at larger distances from 
the star. Figure [4] depicts the situation for a particular 
relative velocity. The avoidance radii of silicate grains 
are about 2 orders of magnitude smaller than those of 
graphite grains. However, the difference in absorption 
coefficients also means that graphite is heated more effi- 
ciently by the radiation field, reaching higher equilibrium 
temperatures and hence contributing more at near-IR 
wavelengths than silicates. Our model shows that for ra- 
diation fields typical of our stars, most of the infrared 
flux, in both MIPS bands, is coming from the carbona- 
ceous grains, with the fluxes produced by silicates being 
at least one order of magnitude smaller, in spite of the 
fact that silicate grains get much closer to the star. 

Figure [4] illustrates another point. The maximum dif- 
ference in avoidance radii between the two compositions 
occurs at the peak of the graphite emission ( a ~ 0.06 
fim). For silicate grains the avoidance radius peaks at a 
slightly larger grain size of about 0.1 fim. As the grain 
size increases to values larger than 0.1 fim, the avoidance 
radii converge to zero. However, our grain size distribu- 
tion, with a sharp cutoff at 0.25 /im, does not have grains 
that are large enough to be accreted onto the stellar sur- 
face. In any case, a limiting avoidance radius is set by 
the sublimation of grains. 

4.4. Total dust mass density 

We have also studied the effect of varying the total 
dust mass density of the cloud, Pcioud- Figure [5] shows 
four different spectral energy distributions for a cloud 
around HD 125162. Again, we have used a grain size of 
0.008 fj,m. and silicate composition. The total gas density, 
ngas, varies between 0.45 cm~'^ and 215 cm~'^, assuming 
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Fig. 5. — Dependence of the dust IR emission on the total dust 
particle density of the cloud. Density values shown correspond to: 
0.45 cm~^ (solid line), 2 cm"-^ (dotted line), 20 cm"-^ (dashed line) 
and 215 cm"-^ (dash-dotted line). 

a gas to dust ratio of 100 and a mean molecular weight of 
1.4 times that of atomic hydrogen. For comparison typ- 
ical diffuse clouds in the galaxy have number densities 
of the order of 10-100 cm~'^. As we go from very diffuse 
clouds to higher densities, keeping all the other param- 
eters constant, the total amount of material increases, 
and similarly, the total flux emitted by the cloud rises. 
The position of the emission peak remains the same, at 
about 55 /zm, since we are not changing the equilibrium 
temperatures of the grains when we increase the density. 

5. DISCUSSION 

5.1. Prevalence of ISM-interaction excesses 

For a given star, our model takes the stellar param- 
eters, the dust optical parameters and the relative ve- 
locity as an input, calculates avoidance radii for grains 
with different sizes using Eq.(2) and finds the total dust 
mass density and size of the cloud that best fit the 
MIPS photometry. This is done for the two limiting ve- 
locities, 15 and 40 km s^^. We perform a non-linear 
least squares fitting to the MIPS photometry. Using 
the reduced y-squared minimisation method proposed by 
iLampton et al.l (Il976[ ). we have derived the parameters 
Pcioud and Tout that best fit the photometry for each of 
our sample stars. Statistical errors for these parameters 
have been estimated at a 90% confidence level. We ob- 
tain number densities for the gas in our clouds, as well as 
total masses of the clouds. The total amount of material 
in the clouds and their gas densities allow a comparison 
with typical values for IS clouds. They are entered in 
Tables 3 and 4, where we also list average metallicities 
for our sample stars. Note that the cloud density values, 
rigas are very low, close to current detection limits. Thus, 
nearly any identified IS cloud can potentially produce an 
excess when an A-star penetrates it. Very diffuse clouds 
with rigas < 1 cm~'^ are more likely to produce an excess 
at higher velocities. 

5.2. Interstellar clouds and debris disks 

Our modelled mid-infrared fluxes are consis tent with 
the Sv itzer photometry of our sample of stars. ISu et al.l 
()2006l ) have shown that the infrared photometry of these 



stars can also be fitted using larger (a few pan) grains dis- 
tributed in debris disks, which means that simple SED 
fitting cannot distinguish between the two phenomena 
(see also iKalas et al.ll2002l ). The uncertainties and vari- 
ations in the ISM-interaction models due to cloud size, 
photometric aperture, and pdust reinforce this point. 

In addition, there is no clear observational distinction 
in the photometric data between the excesses of A Boo- 
tis stars and those of the en tire set of A-stars observed 
with Spitzer (|Su et al.ll2006l ). The incidence of excesses 
is similar, as is the range of 24/70 /im color tempera- 
tures. It is improbable that all A-stars with excesses 
result from interactions with the ISM, but at least in the 
case of 5 Vel^ an ex cess does result from this mechanism 
(|Gaspar et al.|l2008D. Thu s, our m odeling expands on the 
suggestion of iKalas et al.l ()2002[ ) that ISM heating can 
produce far infrared emission that is difficult to distin- 
guish with existing data from a debris disk f or any given 
system. The reflection nebulae discussed in IKalas et al.l 
(|2002.) . however, show no sign of dynamical interaction 
such as bow shocks. The typical mid infrared emission 
(their Table 4) is one to two orders of magnitude larger 
than that from our objects and the derived cloud sizes are 
closer to 100000 AU (their Table 5) as compared to a few 
1000 AU in the interacting systems presented here. We 
thus conclude that the two groups of objects are physi- 
cally different systems. 

Heated ISM would be expected to produce an extended 
distribution of emission by aromatic grains. A more 
definitive test would be to use IFU spectroscopy with 
MIRI on JWST to search for aromatic emission features 
off the star itself (spectra centered on the star will be 
too dominated by the photospheric emission to see such 
features) . 

It is interesting to estimate the probable division be- 
tween debris disk and heated ISM systems. According 
to our results, very diffuse {n = 0.5 — 5 cm"'^) and small 
(''out ~ 0.02 pc) interstellar clouds with total dust masses 
of order 10% of a lunar mass would be enough to produce 
many of the observed infrared excesses. We can make 
a rough estimate of the time an ens emble of stars will 
spend in such clouds from the work of lTalbot fc NewmanI 
(|1977f ). They determined that a typical star like the sun 
has passed through about 135 clouds of 7i(H) > 100 cm~'^ 
and nominal diameter 7 pc in its lifetime, which corre- 
sponds to a probability of being within such a cloud of 
just under 1%. 

The probability for nearby stars to be within such 
a cloud will be lower because of the low density 
of the Local Bubble . Non etheless , it is not zero: 
'Stani mirovic fc HeilesI ()2005[ ) and iBraun fc Kanekarl 
(|2005f ) have detected tiny interstellar clouds within the 
Local Bubble that are closely analogous to our model 
cloud, with diameters of a few thousand AU and densi- 
ties of 5 - 100 cm~'^. It appears that a minority but still 
significant number of infrared excess stars are plausibly 
explained by the ISM-interaction process, given the low 
filling factor of the local ISM. 

Further understanding of this issue requires means to 
identify excesses due to heated ISM dust. One approach 
is to identify bow shock features, which may be charac- 

® and in some cases witho ut clear bow shock morphology such 
as the examples discussed bv IKalas et al.l II2002I) 
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TABLE 3 
Model results for Vj-^i = 15kms~^.^ 



HD No. rigas Tout Mdust ■'"ace M tacc Afacc '■apor,24 ^ rapcr.70 ^avcr 

cm-3 AU Mmooii AU xlO-"MQyr-i yr x10-"Mq AU AU dcx 



11413 


1.00 


2729 


0.091 


16.33 


1.05 


1728 


1.82 


1117 


2617 


-1.07 


30422 


3.12 


1479 


0.045 


10.97 


1.47 


937 


1.38 


859 


2011 




31295 


3.01 


2245 


0.151 


12.97 


1.98 


1422 


2.82 


552 


1294 


-0.91 


110411 


3.71 


1591 


0.065 


12.59 


2.31 


1008 


2.33 


551 


1292 


-0.83 


125162 


5.48 


1402 


0.067 


12.72 


3.47 


888 


3.08 


445 


1043 


-1.35 


183324 


2.80 


1085 


0.011 


13.03 


1.86 


687 


1.28 


881 


2065 


-1.41 



^ The uncertainties in the best-fit parameters (rout and p) are of about 10%, as calculated using the minimisation, 
jigas corresponds to the total number gas density of the cloud, while tacc is the calculated timescale of the star-cloud 
interaction.'' Physical aperture radius at the distance to the star at 24 fim.'^ Physical aperture radius at the distance 
to the star at 70 fim.'^ Average metallicities calculated for the stars based on lHeiterl II2002I ). 



TABLE 4 
Model results for Vj-^i = 40kms~i.'' 



HD No. rtgas Tout Afjust race M tacc Afacc rapcr.24 ^ rapor,70 -^avcr 

cm-3 AU MMoon AU xlO-^^Mgyr-i yr xIO'^^Mq AU AU dex 



11413 


0.30 


4423 


0.117 


2.30 


1.68 


1050 


1.76 


1117 


2617 


-1.07 


30422 


1.24 


2117 


0.053 


1.54 


3.09 


502 


1.56 


859 


2011 




31295 


0.89 


3547 


0.177 


1.82 


3.09 


842 


2.60 


552 


1294 


-0.91 


110411 


1.02 


2425 


0.065 


1.77 


3.33 


576 


1.92 


551 


1292 


-0.83 


125162 


1.25 


2206 


0.060 


1.78 


4.19 


524 


2.19 


445 


1043 


-1.35 


183324 


0.52 


1712 


0.012 


1.83 


1.81 


406 


7.38 


881 


2065 


-1.41 



^ The uncertainties in the best-fit parameters (rout and p) are of about 10%, as calculated using the minimisation, 
rtgas corresponds to the total number gas density of the cloud, while tacc is the calculated timescale of the star-cloud 
interaction.'' Physical aperture radius at the distance to the star at 24 fj.m.'' Physical aperture radius at the distance 
to the star at 70 /im.'^ Average metallicities calculated for the stars based on IHeiteri (2U0% . 



teristic of star-cloud inter actions at the typical values of 
- 30-40 km s-i for v^ei (|Gaspar et all [20081. Careful 
deconvolution of Spitzer imaging (Caspar et al. 2008, in 
preparation) may be capable of finding additional exam- 
ples. Of course, since the avoidance radii derived for our 
program stars are typically larger than the typical debris 
disk outer radius, it is possible that both mechanisms 
could contribute to the excess from a single star. 

Outside the Local Bubble the issue is more difficult, 
both because the number of ISM clouds capable of gen- 
erating excesses is much larger, and because imaging with 
Spitzer is inadequate to identify the characteristic struc- 
ture from an ISM interaction. 

5.3. Contamination of stellar atmospheres 

A star will show A Bootis characteristics only after 
its atmospheric contamination exceeds a threshold. The 
contamination of a stellar atmosphere will be propor- 
tional to the accretio n rate from an interstellar cloud. 
iBondi fc Hovld (|1944f ) solved the equations of motion for 
gas particles that move with uniform speed in the grav- 
itational field of a massive object. They concluded that 
if a star moves across an interstellar cloud with a rela- 
tive velocity Wrei , it accretes material up to an accretion 
radius given by: 



(15) 



M 



cPgas^rel 



(16) 



This radius translates into an accretion rate (the Bondi- 
Hoyle accretion rate) of: 



where pgas is the mass density of gas (analogous to Ugas)- 
This equation assumes that the interstellar gas can be 
treated as a fluid, i.e., the gas kinetic and potential en- 
ergy is efficiently dissipated along the symmetry axis be- 
hind the star. Strictly speaking, this only happens if the 
mean free path is much smaller than the accretion ra- 
dius, a condition that is not met under the conditions in 
our model clouds. The time for the star to move through 
the cloud is short compared to the ionization time scale 
and hence the gas is mostly neutral, so the mean free 
paths are long (although carbon and other metals can be 
largely ionized by the diffuse interste llar UV radiation 
field). However. iBondi fc Hovld (| 19441 ) postulated an in- 
crease in the IS gas density in the wake behind the star, 
since many gas particle trajectories cross on the axis. 
Such behavior is indeed found in multi-dimensio nal nu- 
merica l simulations of Bondi-Hoyle accretion (see lEdead 
(|2004f ) for a summary). The higher density would facil- 
itate accretion. In addition, the best-fit densities from 
our model are " averages" . There is evidence for struc- 
ture in the ISM down to scales o f only t ens of AU (e.g., 
iFrail et al. I (|1994j ). iDeshpandd (|2000f ). iHartquist etall 
(|2003D ). Thus, local very-small-scale density enhance- 
ments may result in accretion at higher rates than im- 
plied by the average densities. 

The details of the accretion are likely, therefore, to 
be far more complex than the Bondi & Hoyle formalism 
applied to the average densities. Nonetheless, we use 
Eqns. (15) and (16) with these densities to make rough 



A Boo Stars 



9 



estimates of accretion rates expected. We estimate the 
accretion radii (race) and the accretion rates (Af) f^or our 
sample stars, assuming they have encountered our model 
cloud. Since we assume that the cloud has a uniform gas 
density, the accretion occurs over the time required for 
the star to cross the cloud. For example, taking a star 
of 2.0 Mq, rigas = 4.0 cm""^, Vrei = 15 km s~^ (from 
Table 3), and rout — 1600 AU (the cloud radius), we 
find an accret i on ra te of 2.3 x 10~^'* Mp, yr~^ (see also 
iGaspar et all ()2008f )). iTurcotte fc Charbonnecii] (|1993l ) 
estimated that an accretion rate > 10~^^ M0 yr~^ will 
yield A Bootis characteristics in a 8000 K main-sequence 
star, which means that our clouds, at least for small ve- 
locities, are in the low density limit of the phenomenon. 

The time for accretion in our small clouds is short; 15 
km s^^ gives a net motion of 3 AU yr~^. Therefore, 
our model cloud is traversed in only 1000 yr and the 
total accreted mass is ~ 2 x 10~^^ M©. Evolutionary 
models show that A-type stars with masses between 2.0 
and 2.5 M© have convection zone masses ranging from 
10^^^ to 10"^" M0. At low velocities, the accreted mass 
is thus comparable to the mass in the convective zone 
and an observable contamination is possible. As we have 
mentioned, two of our stars might have velocities close to 
this lower limit. At larger velocities, close to our upper 
limit of 40 km s~^, the accretion rates drop by two orders 
of magnitude, and for stars with large relative velocities 
it seems unlikely that such a small contamination could 
dominate the abundances in this zone to produce the very 
large anomalies (e.g., two orders of magnitude) seen in A 
Bootis stars. 

Apart from changes in Vj-d, other plausible alternative 
circumstances can produce more favorable results. For 
example, there may be larger clouds, as implied if the 
Local Bubble has a typical self-similar ISM structure. In 
addition, the conditions required by the ISM accretion 
hypothesis become much more plausible outside the Lo- 
cal Bu bble, which in some d irections only extends to ~ 
50 pc (jLallement et al.ll2003| ). 

In summary, we find that the requirements for ade- 
quate accreted matter to produce the A Bootis abun- 
dance signature are significantly more stringent than 
those needed to yield an infrared excess. Alternatively, 
this conclusion can be stated that any star that is actively 
accreting would produce a substantial excess. Nonethe- 
less, this result does not predict that every A Bootis-type 
star should have an infrared excess. The accreted matter 
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from a small IS cloud will remain in the outer layer of the 
star for long enough that the star is likely to have moved 
out of the cloud well before the A Bootis characteristics 
decay away. 

6. CONCLUSION 

We have constructed a geometrical and thermal emis- 
sion model for the interaction of hot stars with diffuse 
interstellar clouds. We find that the predicted flux emit- 
ted by dust in the clouds is consistent with Spitzer pho- 
tometry of A Bootis stars. However, we also find that 
simple SED fitting is not enough to distinguish between 
infrared excesses generated in this manner and those due 
to planetary debris disks. Within the Local Bubble, we 
estimate that the contamination of debris disk samples 
will be low, but further away it may be significant. 

Our model provides order of magnitude estimates for 
the hypothetical cloud parameters. Small {^-^ 0.02 pc) 
and diffuse (~ 0.5 - 5 cm^^) clouds with masses of about 
10% of a lunar mass in solid material are sufficient to 
produce the observed infrared excess around these stars. 
It may be possible to identify examples of such inter- 
actions by observing bow shocks like the one found by 
iGaspar et a l. (2008) for S Velo rum. For these cloud pa- 
rameters, we estimate accretion rates of gas onto the 
stellar surface between 10~^^ and 10~^* Mq yr~^, prob- 
ably too low to produce A Bootis characteristics. How- 
ever, our results suggest that ISM interactions under cer- 
tain favorable conditions (low relative velocity between 
cloud and star, higher than average cloud density, larger 
clouds) may yield both the infrared excesses and abun- 
dance anomalies of such stars. These conditions may be 
relatively common outside the Local Bubble. 
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